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Introduction 

It  is  well  known  that  fractional  differential  equations  both  ordinary  and  fractional  partial  differential 
have  played  an  important  role  in  applications  in  theoretical  modeling  of  scientific  and  engineering  problems. 

The  qualitative  studies  such  as  existence,  uniqueness  and  computation  of  such  dynamical  models  are  useful  in 
applications.  Some  approximate  methods  like  the  radial  functions  method  and  other  methods  are  used  in  [3,  4, 

7,  16]  to  compute  exact  or  approximate  solutions. 

In  solving  nonlinear  problems,  the  monotone  method  combined  with  the  method  of  upper  and  lower 
solutions  is  a  popular  choice,  because  the  existence  of  solutions  by  the  monotone  method  is  both  theoretical  and 
computational.  The  monotone  method  for  various  nonlinear  problems  has  been  developed  in  reference  [8].  The 
monotone  method  (monotone  iterative  technique)  combined  with  the  method  of  lower  and  upper  solutions 
yields  monotone  sequences,  which  converges  to  minimal  and  maximal  solutions  of  the  nonlinear  differential 
equation.  In  many  nonlinear  problems  (nonlinear  dynamic  systems),  the  nonlinear  term  is  the  sum  of  increasing 
and  decreasing  functions.  The  monotone  method  extended  to  such  systems  is  called  the  generalized  monotone 
method.  Here  the  generalized  monotone  method  has  been  extended  to  the  Caputo  fractional  differential  equation 
of  order  q  (where  0  <  q  <  1)  with  an  initial  condition  as  well  as  the  existence  of  coupled  minimal  and  maximal 
solutions  for  such  an  equation  and  a  numerical  example  is  provided  as  an  application  of  the  theoretical  results. 

The  linear  sub  hyperbolic  fractional  partial  differential  equation  in  one  dimensional  space,  that  is  the  qth 
time  order  derivative  is  such  that  1  <  q  <  2,  and  the  linear  super  hyperbolic  fractional  partial  differential 
equation  in  one  dimensional  space,  that  is  the  qth  time  order  derivative  is  such  that  2  <  q  <  3  are  considered  .  In 


the  special  case  when  q=  1,  it  is  the  parabolic  equation  and  when  q  =  2,it  is  the  hyperbolic  equation.  The 
eigenfunction  expansion  method  is  used  to  obtain  representation  forms  for  the  solution  of  the  linear  sub  and  the 
linear  super  hyperbolic  fractional  partial  differential  equations  in  terms  of  the  nonhomogeneous  terms,  initial 
and  boundary  conditions.  See  [9,10,  11,  12,  14]  and  the  references  therein  for  qualitative  study,  applications  and 
numerical  study  of  fractional  differential  equations.  The  representation  forms  obtained  are  useful  in  computing 
numerical  solutions  to  the  sub  and  super  hyperbolic  partial  differential  equations. 


Preliminaries 

In  this  section,  we  recall  some  known  definitions  and  known  results,  which  are  useful  to  develop  our  main 
results.  Initially,  we  recall  some  definitions: 

Definition  2.1.  Caputo  derivative  of  order  q,  when  N  —  1  <  q  <  N  for  t  G  [0,  T], 
is  defined  as 


"DlMt)  =  p^(jo(t  -  sY-^-^u^is)ds). 


In  the  special  case  when  N  =\, 

Definition  2.2.  Caputo  fractional  integral  of  order  q,  when  0  <  q 


< 


r(l-q) 


1,  is  defined  as 
'(s)dsj. 


Definition  2.3.  Riemann  Liouville  fractional  derivative  of  order  q,  when  0<q<lfortE  [0,  T], 
is  defined  as 


D'^  u{f)  = — - —  —  jVt  — s)  ^u(s)ds. 

^  ^  r(i-q)  dt -^0 ^  ^  J 

Definition  2.4.  Riemann  Liouville  right-fractional  integral  of  order  q,  when  0  <  q  <  1  for  t  G  [0,  T], 
is  defined  as 


f-  D  (t  —  s)  ^u(s)ds. 

to  ^  r(q)df>0^  J  ^  J 

Definition  2.5.  The  two  parameter  Mittag-Leffler  function  is  defined  as 


k=0 


r(ak  +  P) 


and  the  single  parameter  Mittag-Leffler  function  is  defined  as  £'Q.(dt“)  = 


Definition  2.6.  let  Q<q<l,  p  —  1  —  q  a  function  0  G  C(J,R)is  a  Cp  continuous  function  if 
(t  —  to)^~‘^0(t)  G  C(Jq,R).  The  set  of  Cp  functions  is  denoted  by  Cp(J,R).  Furthermore,  given  a  function 
0(t)  G  C(Jq,R),  the  function  (t  —  to)^~^0(t)  is  called  the  continuous  extension  of(p(t). 


Consider  the  linear  Caputo  fractional  differential  equation  of  order  q  when  N  —  1  <  q  <  N  with  initial 
conditions  of  the  form 


(  -  lit(t)  =  =  bk  (2.1) 

where  G  M  for  /c  =  0, 1, 2, ... iV  —  1  and  /  G  C[/,  M]  and  J  —  [0,  T].  Also,  denotes  the  kth 

derivative  of  u  with  respect  to  t  at  t  =  0.  If  /(t)  in  (2.1)  is  in  the  space  of  [0,  T] .  The  problem  (2.1)  is 
equivalent  to  the  Volterra  integral  equation 

^  ^  pt  u(s)ds  1  pt  f(s)ds 

f  J  Z.;=0  j,  ppg)  Jo  (t_s)l-q  Y(,q)  -^0  (t-s)!-?' 

The  solution  obtained  is  in  the  space  of  [0,  T],  see  [7]  for  details. 

Consider  the  equation 

(2.5)  (  ^D‘^it)(t)  =  f(t,it(t)  )  +  g{t,u(t)  ),  it(0)  =  ito,  where  0  <  q  <  1, 

f(t,  u{t)  ),  g{t,  u{t)  )  E  Cy  X  R,  R],  f(t,  ii)is  increasing  in  it  on /  and  g{t,  it)  is  decreasing  in  it  on ]. 

Definition.  Let  Vq,  Wq  G  [/,  R] .  Vq,  Wq  are  said  to  be  coupled  lower  and  upper  solutions  of  type  I  of  (2.5)  if 
<  f{t,v{t))  +  g(t,w{t)),  i’(O)  <ito,  and 
>f{t,w{t))  +  g(t,v(t)),  w(0)  >ito. 


Generalized  Monotone  Method  for  Caputo  Fractional  Differential 
Equation  with  an  Initial  Condition 

The  qualitative  study  of  fractional  differential  equations  has  gained  significant  importance  due  to  its 
applications  in  various  branches  of  Science,  Engineering  and  medicine,  see  [1,3,13,14,16,18].  In  [14]  it  has 
been  demonstrated  that  half  order  fractional  differential  equations  give  considerably  better  results  for  certain 
electro  chemical  problems  than  the  classical  approach.  In  addition  [1,13,18]  have  studied  linear  fractional 
reaction  diffusion  equations  which  are  applicable  in  random  walks  and  nanotechnology. 

Consider  the  equation 

(3.1)  (  ^D‘^it)(t)  =  f(t,it(t)  )  +  g{t,u{t)  ),  it(0)  =  ito, 

where  0  <  q  <  l,/(t,it(t)  ),g(t,u(t)  )  G  C[f  x  R,R],  /(t,it)is  increasing  in  it  on /  and  g(,t,u)  is  decreasing 
in  It  on  J. 

Theorem  3.1  Assume  that 

I.  Vq.Wq  E  C^[f,R].  I’o,  Wgare  coupled  lower  and  upper  solutions  of  type  I,  with 

Vq  <  Wgon  J. 

II.  f(t,  u(t)^,  g(t,  u(t)^  E  Cy  X  R,  i?]],  where  /(t,  it)is  increasing  in  it  on /  and  g{t,  it)  is  decreasing  in  it 
on  ]. 

Then  there  exist  monotone  sequences,  i’ji(t)and  Wji(t),  such  that  i’n(t)  p{t)  and  Wji(t)  r(t)  uniformly 
and  monotonically,  where  p(t)and  r(t)  are  coupled  minimal  and  maximal  solutions  of  equation  (3.1)  on  J.  That 
is,  for  any  solution,  u(t),  of  (3.1)  with  i’o(t)  <  ito(t)  <  Wo(t)on  J,  we  get  natural  sequences,  {v^}  and  {w^}, 
satisfying  the  following, 

VQ(t)<V^(t)  <V2(t)  <---<Vn(t)  <it(t)  <  w„(t)  <---<W2(t)  <  Wi(t)  <Wo(t), 


for  each  n  <  1  on  J.  Alsop(t)  <  ii(t)  <  r(t)  on  J. 


The  next  result  uses  intertwined  sequences. 

Theorem  3.2.  Assume  hypotheses,  I.  and  II,  of  the  above  theorem  hold.  Then  for  any  solution,  u(t),  of 
equation  (3.1)  with  r’o(t)  <  Uo(t)  <  Wo(t)on  J,  the  alternating  sequences, 

{v2n>W2n+i}  ^ud  {i^2n+i> Satisfying  intertwined  sequences. 


(3.7)i;o(t)  <  Wi(t)  <•  •  •  <  V2nit)  <  W2n+i(t)  <  u(t)  <  1^2n+l(t)  ^  ^2^(0  <  <  Wo(t), 

for  each  n  >  1  on  /.  This  requires  using  type  (iii)  iterative  schemes, 

CD‘^Vn+i)(t)  =  =  Uoand 

CD‘^Wn+i){t)  =  f(t,Vn{t))  +  p(t,Wn(t)),Wn+i(0)  =  Uq. 

Further,  monotone  sequences,  {i^2n>'^2n+i}  ^^d  {w2n, r’2n+i}»  converge  to  p(t)  and  r(t),  respectively,  on  ]. 
p{t)  and  r(t)  are  coupled  minimal  and  maximal  solutions  of  (3.1),  respectively.  Also,  p(t)  <  ii(t)  <  r(t)  on 

J- 

Theorem  3.4.  If  in  addition  to  hypotheses  I  and  II,  suppose  f(t,  u)  satisfies  left-hand-side  Lipschitz  condition, 
/(t,  112)  —  <  ii(w2  —  u^),U2  >  111,  and  g(t,  u)  satisfies  right-hand-side  Lipschitz  condition, 

p(t,  112)  —  g{t,u^  >  —L2(u2  —  u^),U2  >  111,  where  Vq  <  <  U2  <  Wq  implies  that  p  =  r  =  u,  the 

unique  solution  of  (3.1). 


GENERALIZED  MONOTONE  METHOD  FOR  FRACTIONAL 
REACTION  DIFFUSION  EQUATIONS 


Consider  the  fractional  reaction  diffusion  equation  of  the  form 


a 

d^u  -  =  f(t,x,u)  -f  g(t,x,u)  ,  (t.x)  G 

u(t,0)  =  A(t),  u(t,L)  =  B(t)  (t,x)  G  Tj 
u(0,  x)  —  /i(x)  X  G  n  =  n, 


(4.1) 


where  Cl  —  [0,L],J  =  (0,T],  Q  j  —  J  x  Cl ,  k  >  0  and  Ft’  =  (0, T)  x  dCl. 

Here  is  the  Caputo  Partial  Derivative  with  respect  to  t  of  order  q,  0  <  q  <  1,  which  is  defined  as 


Assume  that/,  g  E  \Qj  x  R,R],  such  that  f(t,  x,  n)  is  non-decreasing  in  n  and  g{t,  x,  n)  is  non-increasing 
in  n,  A(t),5(t)  G  C^[/, i?]  and /i(x)  G  where  0  <  a  <  1. 


In  this  section  we  develop  the  generalized  monotone  method  for  (4.1)  combined  with  coupled  lower  and  upper 
solutions  of  type  I.  The  generalized  monotone  method  yields  monotone  sequences  which  converge  uniformly 
and  monotonically  to  coupled  minimal  and  maximal  solutions  of  (4.1).  Furthermore,  assuming  a  uniqueness 
condition,  we  prove  the  existence  of  the  unique  solution  of  (4.1). 


Theorem  4.1:  Let  (vq,  Wq)  be  coupled  lower  and  upper  solutions  of  (2.1)  such  that  Vq  <  Wq  on  Qj. 
Furthermore,  suppose  f(t,  x,  u)  is  nondecreasing  in  u  on  Qj  and  g{t,  x,  ii)is  nonincreasing  in  Qj  and 

f{t,x,u)—f{t,x,u)  <  L^{u  —  u) 
g{t,x,u)  —  g{t,x,u)  >  —  L2(ii  —  u), 


whenever  u  >  u  on  Qj.  Then  there  exists  monotone  sequences  ^)}  x)}  defined  such  that 

x)  and  Wn(t,  x)  are  solutions  to 


and 


d^v  (t,  x) 

‘'d^Vn(,t,x)  -k — ^ —  =  fit,  x,Vn-i)  +  git, X.W^-f)  ,  (t,x)  G  Qt 

Vnit,0)  =  Ait),  Vnit,L)  =  Bit)  (t,x)  G  Tj 
v^iO,x)  =  /i(x)  X  E  CL  =  CL. 


(4.2) 


d^Wji(t,x) 

d^Wnit,x)  -k - ^ - =  fit,  X,Wn-i)  +  git,  X,Vn-i)  ,  (t,  x)  G  Qj 

w„(t,0)  =  i4(t),  Wnit,L)  =  Bit)  (t,x)  G  T-r 
w„(0,x)  =  /i(x)  X  E  CL  —  CL. 


(4.3) 


respectively.  Also,  these  sequences  converge  uniformly  and  monotonically  to  pit,  x)  and  r(t,  x)  such  that 
ipit,  x),  r(t,  x))  are  coupled  minimal  and  maximal  solutions  of  the  nonlinear  initial  value  problem  (4.1), 
provided  i;o(t, x)  <  ii(t,  x)  <  Wo(t,  x),  where  ii(t,  x)  is  any  solution  of  (2.1).  Furthermore,  p(t, x)  =  r(t,  x)  = 
ii(t,  x),  the  unique  solution  of  (2.1),  on  Qt’. 


A  Representative  Formula  for  the  Linear  Caputo  Fractional  Wave 
Equation 

It  is  well  known  that  fractional  Brownian  motion  has  been  modeled  as  parabolic  stochastic  differential  equation. 
Here  a  representation  form  for  the  solution  of  the  deterministic  one  dimensional  fractional  wave  equation  with 
Caputo  fractional  derivative  of  order  q,  for  1  <  q  <  2  has  been  developed.  For  q  =  1  and  q  =  2,  it 
reduces  to  the  one  dimensional  parabolic  equation  and  one  dimensional  wave  equation,  respectively. 

Consider  the  linear  Caputo  fractional  wave  differential  equation 


^  „  d^u 

dju  -  =  Rit,  x)  ,  it,  x)  E  Qt 

M(t,  0)  =  Ait),  uit,L)  —  Bit)  it,x)  E  Tj 
m(0,  x)  =  hix),UtiO,x)  =  gix)  x  E  Cl, 


(5.1) 


Where  Cl  ^  [0,L],J  =  (0,  T],  Q  ^  =  /  x  Cl ,  >  0  and  =  (0,  T)  x  dCl. 


In  order  to  obtain  a  representation  formula  we  need  to  consider  the  corresponding  homogeneous  equation  with 
homogeneous  boundary  condition.  Using  separation  of  variables  method  we  can  and  the  eigenfunctions 
corresponding  to  the  eigenvalues  The  eigenfunctions  and  the  eigen-  values  vary  depending  on  the  type  of 


boundary  conditions.  In  order  to  get  a  suitable  final  form  of  the  solution  we  normalize  the  eigenfunctions.  For 

(?17r\  ^  TITTX 

—  1  and  the  eigen-  functions  are  (pn(.x)  —  sin—^. 


u(t 


Li  00 

1 

0  in=l 

+Jo  (t  - 

+(_l)n+l. 


Kxo)(pnixo)dxoE^^^i-c^Arit‘')  + 


giXo)(pniXo)dXot‘l 


^Eqqi-c^Ait  -  sy)  Q(is,x)ds +  -  sy  -  sy)A(is)ds 

^  Jo(t  -  -  sy)B(s)ds  . 


Sub  Hyperbolic  Linear  Partial  Fractional  Differential 
Equation  in  One  Dimensional  Space  with  Numerical  Results 


In  this  section  a  representation  form  for  the  sub  hyperbolic  linear  fractional  partial  differential  equation  in  one 
dimensional  space  is  obtained.  It  is  called  sub  hyperbolic  if  the  time  derivative  is  of  order  q  when  1  <  ^  <  2.  It 
is  easy  to  observe  that,  if  <7  =  2,  it  is  a  second  order  linear  hyperbolic  equation  in  one  dimensional  space. 
Consider  the  following  linear  Caputo  fractional  wave  differential  equation: 


d^u 

-  y  =  E{x,  t)  on  Qt, 

u(0,  t)  =  A(t),  u{L,  t)  =  5(t)  on  Ft’, 
it(0,  x)  =  /io(x),  itt(0,x)  =  /ii(x)  on  El, 


(6.1) 


where  12  =  [0,  L],  ]  =  (0,7],  Qt  =  J  x  El,  ci  >  0  and  Fr  =  (0,7)  x  dEl. 


We  develop  a  representation  form  for  (6.1)  which  is  useful  in  computing  the  solution  of  linear  Caputo 
fractional  partial  differential  equations  of  order  q  when  1  <  q  <  2,  using  the  eigenfunction  expansion 

f  Tnix\ 

—  1  and  eigenfunctions  0„(x)  =  sin  (-^1.  This  will  be  useful  in  the 

development  of  the  generalized  monotone  method  to  be  developed  in  future  work.  We  also  assume  that  the 
initial  boundary  condition  satisfy  the  compatability  condition.  Hence,  we  set 

u(t,x)  =  2n=i^n(t)siti^.  (6.2) 

We  must  then  find  h„(t)  for  each  n  >  1.  Taking  the  Caputo  derivative  of  each  side,  Green’s  formula  and  the 
boundary  conditions  yields  the  following  fractional  equation  to  be  solved: 

)  Kit)  = 

With  initial  conditions 

2  y  /nnx\  ,  2  /nnx\ 

hji(O)  —  j  K{x)sin  dx ,  (0)  —  -  J  hi(x)sin  dx. 

Using  Equation  (2.1)  to  solve  (6.3)  substituting  in  for  h„(t)  in  equation  (6.2)  and  interchanging  the  order  of 
integration  and  summation,  we  obtain  the  representation  form: 


1  +  (-l)"S(t)} +  kF(t,x)sm(2^)(ix  (63) 


u(t,  x)  = 


r  00 

II  KXo)(pn(Xo)dXoEq^^{-C 

0  Ln=l 

r  00 

^11  hi(Xo)(pn  (Xo)(Pn(x)dXo 

0  ^n=l 
t  L 


+ 


r  ’Zc^TiTC 

j  ^2  -  5)‘?)i4(s)(p„(x)ds 


r  "Zc^TiK 

+  (-1)^+1  j  (t  -  -  s)'?)5(s)(p„(x)ds. 


Super  Hyperbolic  Linear  Partial  Fractional  Differential 
Equation  in  One  Dimensional  Space  with  Numerical  Results 


Using  the  method  of  successive  approximation  by  setting  u^it)  —  Ylj=Q^  t^  and  taking  the  limit  of 
Uy^{t)  as  m  ^  CO,  the  solution  of  (2.1)  is  given  by 

w(t)  =  +  JgCt  -  sy-^Eci^aiKt  -  sy)fis)ds. 

Next  we  consider  an  eigenvalue  problem  on  one  dimensional  space.  For  that  purpose  consider  the  eigenvalue 
problem 

'll"  +  Xll  —  0 

with  mixed  homogeneous  boundary  conditions  of  the  form: 

5y(a)  =  cciu(i)  +  —  O,fori  =  0,1. 

The  value  of  X  for  which  the  nontrivial  solution  for  the  boundary  value  problem  exists  is  called  the  eigenvalue 
denoted  as  Xj^  and  the  corresponding  nontrivial  solution  is  called  the  eigenfunction  denoted  as  (Pnipd)  ■  The  main 
result  is  developed  using  the  eigenvalue  problem  with  Direchlet  boundary  conditions, 

m(0)  =  u{L)  =  0. 


In  this  section  a  representation  form  for  the  super  hyperbolic  linear  fractional  partial  differential  equation  in  one 
dimensional  space  is  obtained.  It  is  called  super  hyperbolic  if  the  time  derivative  is  of  order  q  when  2  <  <5^  <  3.  It 
is  easy  to  observe  that,  if  <7  =  2,  it  is  a  second  order  linear  hyperbolic  equation  in  one  dimensional  space. 
Consider  the  following  linear  Caputo  fractional  wave  differential  equation: 


^  „  „ d^u 

d^u  -  =  F(x,  t)  on  Qt, 

ii(0,  t)  =  A(t),  u{L,  t)  =  5(t)  on  Ft^, 
ii(0,x)  = /io(x),  iit(0,x)  =  /ii(x),  iitt(0,x)  =  /i2(x)  onfi, 


(7.1) 


where  12  =  [0,  L],  ]  =  (0,7],  Qt  =  J  x  Cl,  cz  >  0  and  Fr  =  (0,7)  x  dfi.  We  also  assume  that 
7l(t),5(t)  G  C3[/,F],  and F(t,x)  G  C^'^[Qr,R]. 


Assume  the  following  compatability  conditions  relative  to  the  initial  and  boundary  conditions: 


Fo(0)  =  A(0)  =  0 


Fo(L)  =  5(0)  =  0 


/Ii(0)  =^'(0)  =  0, 

/i2(0)=^”(0)  =  0, 


/ii(L)  =  5''(0)  =  0 
/i2(L)  =5”(0)  =  0. 


This  is  needed  when  differentiating  the  integral  representation  of  the  terms  related  to  the  boundary 
conditions  and  when  the  problem  is  reduced  to  homogeneous  boundary  conditions.  As  in  section  3,  we  set 
u(t,  x)  =  As  before,  the  expression  for  by^(t)  is  found  for  each  n  >  1  by  using  Lagrange’s 

identity  along  with  Green’s  formula  and  the  boundary  conditions  and  substitute  back  into  the  expression  for 

ii(t,  x).  Changing  the  order  of  integration  and  summation,  we  obtain  the  solution  : 

u(t,x)  = 


00  2 

hjixo)  (pn  (xo)dxot^'5qj+i(-c^ 

n=l j=0 
L 


+  ^j  Ari(.t  -  sy)[F(,S,Xo)(PniXo)(Pn(x)dXo]ds 

0  0 

t  00 

r  "Zc^Tni 

+  J  ^  ^2  -  sY)A{s)(Pn{x)ds 

0  n=l 


+  (-l) 


n+1 


L 

/ 


2c^nn 


(t  -  -  s)'?)5(s)(p„(x)ds. 


0 


This  representation  formula  for  the  solution  of  (7.1)  provides  a  methodology  for  computing  the  solution  for 
different  nonhomogeneous  terms  and  different  initial  and  boundary  conditions.  The  convergence  of  the  infinite 
series  involves  the  Mittag-Leffler  functions  is  yet  to  be  established  even  when  0  <  q  <  1.  However,  if  the 
boundary  conditions  are  homogeneous  and  the  initial  and  nonhomogeneous  terms  are  chosen  such  that  there  are 
a  finite  number  of  Mittag-Leffler  functions,  then  the  solution  can  be  numerically  computed  using  MAPLE. 


Riemann  Liouville  and  Caputo  Fractional  Differential  and 
Integral  Inequalities 

In  this  section,  results  for  coupled  fractional  and  ordinary  integral  inequalities  where  the  nonlinear  function  is 
the  sum  of  an  increasing  and  a  decreasing  function  have  been  developed.  Also,  the  corresponding  coupled 
fractional  and  ordinary  differential  inequality  results  have  been  developed  without  requiring  the  increasing  or 
decreasing  nature  of  the  nonlinear  function.  In  fact,  results  have  been  developed  for  coupled  differential  and 
integral  inequalities  for  both  ordinary  and  fractional  equations.  These  results  are  useful  in  the  qualitative  study 
of  ordinary  and  fractional  dynamic  systems  of  both  Riemann  Liouville  and  Caputo  forms. 

For  the  sake  of  simplicity,  the  obtained  results  will  be  on  the  interval  ]  =  (to,  T].  Furthermore, /o  =  (to,  T]  =  /. 

Theorem  8.1.  Let  i;,  w  G  Cp\j  ,R\,f ,  g,E  C[/o  x  R,R\  and  satisfy  the  following  coupled  integral  inequalities 
(i)v(t)  <  +  ^[f{s.v(s))  +  5'(s,w(s)))ds, 


(n)w(t)  >  +  (f(s,w(s))  +  g(s,v(s)))ds, 

(in)  f(t,  u)  is  nondecreasing  in  ufor  each  t  G  Jo,  and  f(t,  u)  satisfies  the  one  sided  Lipschitz  condition  of  the 
form 

f(t,u^)  -  f{t,U2)  <  L(u^-  U2), 
for  some  L  >  0,  whenever  >  112 ; 

(iv)  g(t,  u)  is  non  increasing  in  ufor  each  t  G  Jo  D  and  g(t,  u)  satisfies  the  one  sided  Lipschitz  condition  of 
the  form 

g(t,u^)  -  g(t,U2)  >  -M(u^  -  U2), 
for  some  M  >  0,  whenever  >  112.  Then  Uq  <  Wq  implies  v(t)  <  w(t)  on  J. 

Consider  the  Riemann-Liouville  fractional  differential  equation  of  the  form: 

D^u(t)  =  f(t,u(tf)  +  g(t,u(t)),  r(q)(t  -  (8.1) 

where  f,g  G  C[/o  xR,R]. 

The  next  result  is  the  differential  inequality  result  relative  to  the  Riemann-Liouville  fractional  differential 
equation  (8.1). 

Theorem  8.2.  Let  r;,  w  G  Cp\J,R],f,g,E  C\Jo  x  R,RJ  and  satisfy  the  following  coupled  differential 
inequalities 

(i) D^{v(t))  <  f{t,v(t))  +  g{t,w(t)),  h(q)(t  -  toy-‘iv(t)\t=t^  < 

(ii) D^{w(t))  >  f{t,w(t))  +  g{t,v(t)),  r(q)(t  -  to)^"‘?w(t)|t^t^  > 

(Hi)  f(t,u)  satisfies  the  one  sided  Lipschitz  condition  of  the  form 

-  f(t,U2)  <  L(u-^-  U2), 

for  some  L  >  0,  whenever  u-^  >  U2', 

(iv)  g(t,u)  satisfies  the  one  sided  Lipschitz  condition  of  the  form 

g(t,u^)  -  g(t,U2)  >  -M(u^  -  U2), 
for  some  M  >  0,  whenever  u-^  >  U2.  Then  Vq  <  Wq  implies  v(t)  <  w(t)  on  J. 

Here  the  Caputo  fractional  differential  and  integral  inequalities  are  developed. 

Consider  the  Caputo  fractional  differential  equation  of  order  q  where  0  <  q  <  1,  of  the  form  : 

^D^u(f)  —  A  u(t)  +  f(f),  u(to)  —  uo,  (8.2) 

where  f,g  E  C\Jo  x  R,R]. 

The  integral  representation  of  (8.2)  is  given  by  equation 

u(t)  =  uo  +  Lq,i(l(t  -  to)‘?)  +  -  5)^1-^  -  to)‘?)/(s)ds,  (8.3) 

where  L(q)  is  the  Gamma  function.  Note,  the  above  holds  true  if  in  place  of  equality  we  have  less  than  or  equal. 


Now  consider  the  Caputo  fractional  differential  equation 


(8.4) 


^D^u{t)  =  f{t,u)  +  g{t,u),  ii(to)  =  no, 

Where/, 5-  G  C[/o  x  R,R\. 

Definition  8.1.  The  functions  v,w  E  C^([to,  T],  R)  are  called  natural  lower  and  upper  solutions  of  (8.4)  if: 

<  f(t,v)  +  g(t,v),  v(to)  <  no, 

and 

^D‘^w(t)  >/(t,w)+  g(t,w),  w(to)  >  no. 

Definition  8.2.  n,w  G  C^([tQ,T],R)  are  called  coupled  lower  and  upper  solutions  of  type  I  of 

(8.4)  if: 


^D^v(f)  </(t,  n)+  g(t,w),  n(to)  <  no, 

and 

^D^w(t)  >/(t,w)+  g(t,v),  w(to)  >  no. 

One  can  easily  prove  the  existence  of  the  solution  of  (8.4)  on  the  interval  [to,  rjwhen  natural  lower  and  upper 
solution  for  (8.4)  are  known  for  which  v(t)  <  w(t).  Coupled  minimal  and  maximal  solutions  of  (8.4)  on  the 
interval  [tg,  T]  can  be  computed  without  any  extra  assumptions  using  the  generalized  monotone  method  when 
coupled  lower  and  upper  solutions  of  type  I  of  (8.4)  are  known  with  v(t)  <  w(t). 

Next  we  have  the  differential  inequality  result  for  coupled  lower  and  upper  solutions  of  type  I: 

Theorem  8.3.  Let  v,w  E  C^(J,  R),  Where  f,g  G  C\J  x  R,R]  are  coupled  lower  and  upper  solutions  of  type 
I,  such  that 

(i)  f(t,  u)  satisfies  the  one  sided  Lipschitz  condition  of  the  form 

f(t,u^)  -  f{t,U2)  <  L(u^-  U2), 

for  some  L  >  0,  whenever  >  112 ; 

(ii)  g(t,  u)  satisfies  the  one  sided  Lipschitz  condition  of  the  form 

g(t,u-f}  -  g(t,U2)  >  -  U2), 

for  some  M  >  0,  whenever  u-i  >  U2.  Then  vfto)  <  w (Iq)  implies  that  v(t)  <  w(t),onJ. 

Remark  5.1.  Caputo  integral  inequalities  results  similar  to  Theorem  8.1  can  also  be  developed.  The  integral 
inequality  results  need  an  extra  assumption  which  is  not  required  for  the  corresponding  scalar  differential 
inequality  results.  Thus,  the  differential  inequality  results  are  useful  in  iterative  methods  like  the  monotone 
quasilinearization  methods.  They  are  useful  for  proving  uniqueness,  continuous  dependency  on  the  initial 
condition  and  also  in  determining  the  order  of  convergence. 

Numerical  Graphs 

Consider  a  special  case  of  the  Logistic  model  of  the  form 
(4.1)  =  0.99u(t)  -  u2(t),u(0)  =  i. 
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Figure  2.  Dashed:  {r„}.  Solid: 


Conclusion 


We  have  developed  generalized  monotone  method  for  nonlinear  fractional  differential  systems  and  fractional 
reaction  diffusion  equations  via  coupled  lower  and  upper  solutions  of  type  1 .  The  advantage  of  the  generalized 
monotone  method  is  that  each  component  of  the  iterates  are  scalar  fractional  differential  equation  whose 
solutions  are  easy  to  compute  compared  to  computing  the  solution  corresponding  linear  system. 

Also,  This  representation  formulas  provide  a  methodology  for  computing  solutions  for  different 
nonhomogeneous  terms  and  different  initial  and  boundary  conditions. 
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